Strong exciton-plasmon coupling in semiconducting carbon nanotubes 
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We study theoretically the interactions of excitonic states with surface electromagnetic modes 
of small-diameter ( < 1 nm) semiconducting single-walled carbon nanotubes. We show that these 
interactions can result in strong exciton-surface-plasmon coupling. The exciton absorption line- 
shapes exhibit the line (Rabi) splitting ~ 0.1 — 0.3 eV as the exciton energy is tuned to the nearest 
interband surface plasmon resonance of the nanotube. We expect this effect to open a path to new 
optoelectronic device applications of semiconducting carbon nanotubes. 
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Single-walled carbon nanotubes (CNs) are quasi-one- 
dimensional (ID) cylindrical wires consisting of graphcnc 
sheets rolled- up into cylinders with diameters ~ 1 — 10 nm 
and lengths ~ 1 — 10 4 ^m 0, 0, CNs are shown 
to be useful for miniaturized electronic, electromechani- 
cal, chemical and scanning probe devices and as materi- 
als for macroscopic composites [3] • The area of their po- 
tential applications has been recently expanded towards 
nanophotonics and optoelectronics p, @| after the exper- 
imental demonstration of controllable single-atom incap- 
sulation into single- walled CNs Q. 

For pristine (undoped) single-walled CNs, the numeri- 
cal calculations predicting large exciton binding energies 
(~ 0.3 — 0.6 eV) in semiconducting CNs [1] and even in 
some small-diameter (~ 0.5 nm) metallic CNs fol- 
lowed by the results of various measurements of the exci- 
tonic photoluminescence [l(J 11. 12|. have become avail- 
able. These works, together with other reports investi- 
gating the role of effects such as intrinsic defects [Tll j . 
exciton-phonon interactions I12L Il3l . Il4|. external mag- 
netic and electric fields [IH, Il6| . reveal the variety and 
complexity of the intrinsic optical properties of carbon 
nanotubes. 

Here we develop a theory for the interactions between 
excitonic states and surface electromagnetic modes in 
small-diameter (< 1 nm) semiconducting single-walled 
CNs. The approach is based on our recently developed 
Green function formalism to quantize an electromag- 
netic field in the presence of quasi-lD absorbing bod- 
13| . We show that such interactions can re- 
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suit in a strong exciton-surface-plasmon coupling due to 
the presence of the low-energy (~ 0.5 — 2 eV) weakly- 
dispersive interband plasmon modes [2Cj| and the large 
exciton excitation energies ~ 1 eV [2lj in small-diameter 
CNs. This may be used for new CN optoelectronic device 
applications. 

We consider the vacuum-type electromagnetic interac- 
tion of an exciton with surface electromagnetic fluctua- 
tions of a single- walled semiconducting carbon nanotube. 



No external electromagnetic field is assumed to be ap- 
plied. Since the problem has the cylindrical symmetry, 
the orthonormal cylindrical basis {e r ,e v ,e z } is used with 
the vector e z directed along the nanotube axis. The total 
Hamiltonian of the coupled exciton-photon system is of 
the form (we use Gaussian units) 



H = H f 



Hi, 



(1) 



where the three terms represent the free field, the free ex- 
citon, and their interaction, respectively. More exp licitly, 
the second quantized field Hamiltonian is [Tfl, ll2| 
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Hf = V / dujhuj p{n,u)f(n,uj) , 
- Jo 
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where the scalar bosonic field operators P(n, w) and 
/(n, cj) create and annihilate, respectively, the surface 
electromagnetic excitation of frequency w at an arbi- 
trary point n = R n = {Rcn , <Pn, z n\ associated with 
a carbon atom (representing a lattice site) on the sur- 
face of the CN of radius Rcn- The summation is made 
over all the carbon atoms, and in the following it is re- 
placed by the integration over the entire nanotube sur- 
face according to the rule J3 n - • • = (l/So)fdR n . . . = 

(l/S )f^d<p n R C N JZ* dz "- ■ ■ > where S ° = (3V3/4)& 2 is 
the area of an elementary equilateral triangle selected 
around each carbon atom in a way to cover the entire 
surface of the nanotube with 6= 1.42 A being the carbon- 
carbon interatomic distance. 

The second quantized Hamiltonian of the free exciton 
is of the form (see, e.g., Ref. [22j) 



H F =J2Ef(n)Bl +mJ B mJ = 



(3) 

where the operators B^ a , and B n j create and annihi- 
late, respectively, an exciton with energy Ef(n) at the 
lattice site n of the CN surface. The index / (^ 0) refers 
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to the internal degrees of freedom of the exciton. Alterna- 
tively, B^j = En^l ye 4k n / y/N creates (N is the number 

of the lattice sites) and B^j = annihilates the 

exciton with quasi-momentum k = {k Vl k z } and energy 

£/( k ) = E„£/(n)e- lk - n /iV = ^£ + h 2 k 2 /2M ex , where 



E 



if). 
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exc — is the /-internal-state exciton excitation 

energy, Ejf is the corresponding binding energy of the 
exciton, M ex ~m e + mh is its effective translational mass, 
and Eg is the nanotubc band gap. The two equivalent 
free-exciton Hamiltonian representations are related to 
one another via the obvious orthogonality relationships 

£ n£r i(k-k').n /iV = ^ £ ke -i(n- m) .k /7V = 

with the k-summation running over the first Brillouin 
zone of the nanotube. The bosonic field operators in 
Eq. ((2]) are transformed to the k-representation in the 
same way. 

The most general (non-relativistic, electric dipole) 
exciton-photon interaction on the nanotube surface can 
be written in the form (see Refs. [181 Il9j) 



H 



(i) 



H 



(2) 
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A(n)l +^d n - V n £(n). 



where p n = Ylf(Q\Pn\f)B n j + h.c. is the total electron 
momentum operator at the lattice site n under the opti- 
cal dipole transition resulting in the exciton formation at 
the same site, d n = ^2f(0\d a \f)B a j + h.c. is the corre- 
sponding dipole moment operator [related to p n via the 
equation (0|p n |/) = im e Ef(n)(0\d n \f)/he]. The vector 
potential operator A(n) (the Coulomb gauge is assumed) 
and the scalar potential operator (p(n) represent, respec- 
tively, the nanotube's transversely polarized surface elec- 
tromagnetic modes and longitudinally polarized surface 
electromagnetic modes which the exciton interacts with. 
We express them in terms of our earlier developed elec- 
tromagnetic field quantization formalism in the presence 
of quasi-lD absorbing bodies to obtain IE n| 



A(n) = / du — E (n, w) + h.c. = > / dtu (5) 
Jo ^ ^ Jo 

x-\/nfiu)Re<Jzz{RcN,u) ^G zz (n, m, w)/(m, u) + h.c. 



and 
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Vn^(n) = / dwE (n,w) + h.c. = V / du> (6) 
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—5- nhujRea zz (R C n , u>) "G zz (n, m, w)/(m, w) + h.c. 
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with the total electric field operator given by E(n) = 

f™du\E ± (n,w)+ti(n,w)]+h.c., x CU)G„(n,m,w) being 
the zz-component of the transverse (longitudinal) Green 



tensor (with respect to the first variable) of the electro- 
magnetic subsystem, and o zz (Rcn representing the 
CN dynamic surface axial conductivity per unit length. 

Equations (P)-© form the complete set of equations 
describing the exciton-photon coupled system on the CN 
surface in terms of the electromagnetic field Green tensor 
and the CN surface axial conductivity. The conductivity 
is found beforehand from the realistic band structure of 
a particular CN. The Green tensor is derived by expand- 
ing the solution of the Green equation in cylindrical co- 
ordinates and determining the Wronskian normalization 
constant from the appropriately chosen boundary condi- 
tions on the CN surface (see Refs. [H EE EE 0| ) . More 
details on these calculations will be given elsewhere. 

It is important to realize that the transversely polar- 
ized surface electromagnetic mode contribution to the 
interaction Hamiltonian from Eq. ([4]) (first term) is 
negligible compared to the longitudinally polarized sur- 
face electromagnetic mode contribution (second term). 
The point is that, because of the nanotube quasi-one- 
dimensionality, the exciton quasi-momentum vector and 
all the relevant vectorial matrix elements of the momen- 
tum and dipole moment operators are directed predomi- 
nantly along the CN axis (the longitudinal exciton) . This 
prevents the exciton from the electric dipole coupling to 
transversely polarized surface electromagnetic modes as 
they propagate predominantly along the CN axis with 
their electric vectors orthogonal to the propagation di- 
rection. The longitudinally polarized surface electromag- 
netic modes are generated by the electronic Coulomb 
potential (see, e.g., Ref. (HI), and therefore represent 
CN surface plasmon excitations. These have their electric 
vectors directed along the propagation direction. They 
do couple to the longitudinal excitons on the CN sur- 
face. Such modes were observed in Ref. [23]. They occur 
in CNs both at high energies (well-known 7r-plasmon at 
~6 eV) and at comparatively low energies of ~0.5— 2 eV 
[plasmon modes associated with (transversely quantized) 
interband electronic transitions]. 

To obtain the dispersion equation of the coupled 
exciton-plasmon excitations, we utilize Bogoliubov's ca- 
nonical transformation technique (see, e.g., Ref. [HI]) and 
diagonalize the Hamiltonian H])-© exactly. This results 



H = 



^(k)|t (k)4(k)+£o 



(7) 



k, fi=l,2 



where |„(k) = £/K(lc, Uf)B kJ - » M (k,w/)fll k)/ ] + 
J dwfu^k, w)/(k, w) — u„(k, u))p(— k, u)] annihilates 
and £t(k) = [^(k)]* creates the coupled exciton-plasmon 
excitations. Here /x enumerates the two exciton-plasmon 
branches, u M and are the (appropriately chosen) 
canonical transformation coefficients, u>f = Ef/h. The 
energy Eq is the "vacuum" energy with no exciton- 
plasmons excited in the system, and ftw M (k) is the 
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exciton-plasmon energy given by the solution of the fol- 
lowing (dimensionless) dispersion equation 
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£ f rfi \ g(g) 
/ dxxT J (x) — 
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(8) 



Here x = huj/2"fo, %n = ^^(k)/27o, £/ = i?/(k)/27o with 
7o = 2.7 eV being the carbon nearest neighbor overlap 
integral entering the CN surface axial conductivity a zz . 
The function t f (x) = {2~f /h) 2 (4\d f \ 2 x 3 /3hc 3 ), where 
df = 5Z n (0|d n |/), represents the (dimensionless) spon- 
taneous decay rate, and 



p(x) = 



35*0 



Re- 
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AiraR 2 CN "~ a zz (x) 



(9) 



stands for the surface plasmon density of states (DOS) 
responsible for the decay rate variation due to the cou- 
pling to plasmon modes. Here a = e 2 /Tic = 1/137 is 
the fine-structure constant and a zz — 2itTi<T zz / e 2 is the 
dimensionless CN surface axial conductivity per unit 
length. The conductivity factor equals Re(l/cr zz ) = 
— (8n /ujRc n )Im[l / (e zz — 1 ) in view of the Drude relation- 
ship a zz = —iuj(e zz — l) /AirSpr, where e zz is the longitudi- 
nal (along the CN axis) dielectric function, S and pr are 
the surface area of the tubule and the number of tubules 
per unit volume, respectively 17, [HI, EH]. This relates 
very closely the surface plasmon DOS function ^ to 
the loss function — Im(l/e) measured in Electron Energy 
Loss Spectroscopy (EELS) experiments to determine the 
properties of collective electronic excitations (ioj . 

The calculated plasmon DOS functions are shown in 
the left panels of Fig. Q] (a), (b) and (c), along with 
the functions Kea zz (x) and lma zz (x), for the (11,0), 
(10,0) and (8,0) CN we study here (R CN =0.43, 0.39 
and 0.31 nm, respectively). In all graphs the lower di- 
mensionless energy limits are set up to equal the lowest 
bright exciton excitation energy [E exc =1.21, 1.00 and 
1.55 eV for (11,0), (10,0) and (8,0) CN 0, respectively]. 
Peaks of the plasmon DOS are seen to coincide in energy 
with the zeros of the imaginary conductivities (or the 
zeros of the real dielectric functions), clearly indicating 
the plasmonic nature of the CN surface excitations un- 
der consideration. Peaks broaden as the CN diameters 
decrease. This is consistent with the stronger hybridiza- 
tion effects in smaller-diameter CNs [2(| . The CN dielec- 
tric response is calculated in the following way. First, we 
adapt the nearest-neighbor non-orthogonal tight-binding 
approach [27| to determine the realistic band structure 
of each CN. Then, the room-temperature longitudinal di- 
electric function and conductivity are calculated within 
the random-phase approximation with electronic dissipa- 
tion processes included in the relaxation-time approxima- 
tion (electron scattering length of 130i?cw was used [I3|)- 

We further take advantage of the peak structure of 
p(x) and solve Eq. ((5|) analytically using the Lorentzian 
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FIG. 1: (Color online) (a),(b),(c) Surface plasmon DOS and 
conductivities (left panels), and lowest bright exciton disper- 
sion when coupled to plasmons (right panels) in (11,0), (10,0) 
and (8,0) CN, respectively. Dimensionless energy is defined as 
[Energy] /270; see text for dimensionless quasimomentum. 



Ax 2 p ] 



approximation with p{x) Ri p(x p )Ax 2 /[(x — x p ) 2 
where x p and Ax p are, respectively, the position and the 
half-width-at-half-maximum of the plasmon resonance 
(Fig. [TJ left panels) closest to the lowest bright exciton 
excitation energy in the same nanotube. We obtain 



X1.2 




l^(e 2 f -x 2 p ) 2 + F p e f , 



(10) 



F p = F(x p )Axp(l- Ax p /nx p ), F{x p ) = 2x p T{ j (x p )p{x p ). 
The dispersion curves thus obtained are shown in the 
right panels of Fig. Q] (a), (b) and (c) as functions of the 
longitudinal dimensionless quasi-momentum. All graphs 
demonstrate a clear anticrossing behavior of the two 
exciton-plasmon branches with the (Rabi) energy split- 
ting ~ 0.1 eV. In these calculations, we estimated the 
interband transition matrix element in Tq(x p ) from the 
relationship |d/| 2 = 3h\ 3 /4T™ d f28j, where T r e f is the 
exciton intrinsic radiative lifetime, A = 2irch/E with 
E = E exc + (2irh/3b) 2 t 2 /2M ex being the total energy 
and — 1 < t < 1 representing the longitudinal dimension- 
less quasi-momentum of the exciton. We used the lowest 
bright exciton parameters reported in Ref. 
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We also derive the (dimensionless) exciton absorption 
lineshape function I(x) in the vicinity of the plasmon 
resonance in the way this was done in Ref. [5] for the 
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FIG. 2: (Color online) (a),(b),(c) Exciton absorption line- 
shapes as the exciton energies are tuned to the nearest plas- 
mon resonance energies (vertical dashed lines here; see Fig.[T] 
left panels) in the (11,0), (10,0) and (8,0) CN, respectively. 

optical absorption by atomically doped CNs. In doing 
so, we take into account the exciton-phonon scattering 
in the relaxation time approximation. We obtain 



I(x)=I (e f ) 



{x-e f f + A\ 



+ (x-e f )*A 



(11) 



where I (e f ) = r£(g/ )p(e/)/27r, X f = ^AnAx p Tp jej) , 
Ai = yj2Ae f + Ax p , and A 2 = y/{As f ) 2 + A\ with 
Ae f — h/2jQT p h being the exciton energy broadening due 
to the phonon scattering with the relaxation time T p h- 
The calculated exciton absorption lineshapes for the CNs 
under consideration are shown in Fig. [3] (a) , (b) and (c) 
as the exciton energies are tuned to the nearest plasmon 
resonances. We used r p h = 30 fs as reported in Ref. 
The clear line (Rabi) splitting effect is seen to be of 
the order of 0.1 — 0.3 eV, indicating the strong exciton- 
plasmon coupling with the formation of the mixed sur- 
face plasmon-exciton excitations. The splitting is larger 
in the smaller diameter nanotubes, and is not masked by 
the exciton-phonon scattering. 

In summary, we have shown the strong exciton-surface- 
plasmon coupling effect with characteristic Rabi splitting 
^0.1—0.3 eV in small-diameter (<1 nm) semiconducting 
CNs. This is almost as large as typical exciton binding 
energies in such CNs (^0.3—0.6 eV Q), and of the same 
order of magnitude as the exciton-plasmon Rabi splitting 
in organic semiconductors (~ 180 meV [1^]). Also, this is 
much larger than the exciton-polariton Rabi splitting in 
semiconductor microcavities (~ 140—400 /ieV [30j), or the 
exciton-plasmon Rabi splitting in hybrid semiconductor- 
metal nanoparticle molecules [3l|. However, the forma- 
tion of the strongly coupled exciton-plasmon states is 
only possible if the exciton total energy is in resonance 
with the energy of an interband surface plasmon mode. 



That is not always the case experimentally. To realize the 
strong coupling effect, the exciton energy should be tuned 
to the nearest plasmon resonance in ways used for the 
excitons in semiconductor microcavities — thermally (3pj 
(by elevating sample temperature), and/or electrostati- 
cally [13] (via the quantum confined Stark effect with an 
external electrostatic field applied perpendicular to the 
nanotube axis). This opens paths to new optoelectronic 
device applications of semiconducting CNs. 
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